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1 Introduction 



Reaction-diffusion systems, is a well-studied area. People have studied reaction- 
diffusion systems, using analytical techniques, approximation methods, and sim- 
ulation. The approximation methods may be different in different dimensions, 
as for example the mean field techniques, working good for high dimensions, 
generally do not give correct results for low dimensional systems. A large frac- 
tion of analytical studies, belong to low-dimensional (specially one-dimensional) 
systems, as solving low-dimensional systems should in principle be easier. [1—13] . 

The term exactly-solvable have been used with different meanings. In [14- 
16], intcgrability means that the iV-particlc conditional probabilities' S-matrix is 
factorized into a product of 2-particle S-matrices. In [17-26], solvability means 
closedness of the evolution equation of the empty intervals (or their generaliza- 
tion. 

In [27], a ten-parameter family of reaction-diffusion processes was intro- 
duced for which the evolution equation of n-point functions contains only n- or 
less- point functions. We call such systems autonomous. The average particle- 
number in each site has been obtained exactly for these models. In [28,29], this 
has been generalized to multi-species systems and more-than-two-site interac- 
tions. 

Among the important aspects of reaction-diffusion systems, is the phase 
structure of the system. The static phase structure concerns with the time- 
independent profiles of the system, while the dynamical phase structure concerns 
with the evolution of the system, specially its relaxation behavior. In [30- 
33], the phase structure of some classes of single- or multiple-species reaction- 
diffusion systems have been investigated. These investigations were bases on 
the one-point functions of the systems. 

Here we want to study the two-point functions of autonomous single-species 
one-dimensional reaction-diffusion systems. Throughout this study, the initial 
condition of the system is taken to be translationally- invariant, so that it re- 
mains translational-invariant during the evolution. The two-point function for 
such systems is obtained, and it is shown that it exhibits a non-trivial dy- 
namical phase structure. In section 2, the evolution equation of the two-point 
function is obtained. In section 3, this equation is solved and the corresponding 
energy-spectrum is obtained. In section 4, the parameter space of the system 
is analyzed. In section 5, the dynamical phase structure of the system (the 
different phase regions in the parameter space) is investigated. Finally, section 
6 is devoted to a one-parameter example family, which can be in all five phases. 

2 Evolution equations of the one- and two-point 
functions 

Consider a one-dimensional periodic lattice, every point of which is empty or 
contains one particle. Let the lattice have L + 1 sites. The observables of such 
a system are the operators N", where i with 1 < i < L + 1 denotes the site 
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number, and a — 0, 1 denotes the hole or the particle: N® is the hole (vacancy) 
number operator at site i, and is the particle number operator at site i. One 
has obviously the constraint 

s a N? = 1, (1) 

where s is a covector the components of which (s Q 's) are all equal to one. The 
constraint Q, simply says that every site is either occupied by one particle or 
empty. A representation for these observables is 

2Vf := !®- - -®l ®N a ® 1 <8> • • • ® 1 , (2) 

i— 1 L+l-i 

where N a is a diagonal 2x2 matrix the only nonzero element of which is the 
a'th diagonal element, and the operators 1 in the above expression are also 2x2 
matrices. It is seen that the constraint JU can be written as 

s-N=l, (3) 

where N is a vector the components of which are N°"s. The state of the system 
is characterized by a vector 

P G V(8)---<8>V, (4) 
. ' 

L + l 

where V is a 2-dimensional vector space. All the elements of the vector P are 
nonncgative, and 

S-P = l. (5) 

Here S is the tensor-product of L + 1 co vectors s. 

As the values of the number operators N" 1 are zero or one (and hence iV"'s 
are idempotent), the most general observable of such a system is the product 
of some of these number operators, or a sum of such terms. Moreover, the 
constraint Q shows that the two components of Nj are not independent, so, 
one can express any function of Nj in terms of 

ni := a • Nj, (6) 

where a is an arbitrary covector not parallel to s. Our aim is to study the 
evolution of the two-point functions constructed by n^s. 
The evolution of the state of the system is given by 

P = HP, (7) 

where the Hamiltonian 7i is stochastic, by which it is meant that its nondiagonal 
elements are nonnegative and 

S7i = 0. (8) 
The interaction is nearest-neighbor, if the Hamiltonian is of the form 

L + l 

H = J2 H i,i+i, ( 9 ) 

i=l 
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where 

ff. M+ i := 1 ® ■ ■ ■ ® 1 (gig g> 1 g> ■ ■ ■ g> 1 . (10) 

i— 1 L— « 

(It has been assumed that the sites of the system are identical, that is, the 
system is translation-invariant. Otherwise H in the right-hand side of IjlUfl 
would depend on i.) The two-site Hamiltonian H is stochastic, that is, its non- 
diagonal elements are nonnegative, and the sum of the elements of each of its 
columns vanishes: 

(s(gs)# = 0. (11) 

Here H is a 4 x 4 matrix (as the system under consideration has two possible 
states in each site and the interactions are nearest neighbor). The non-diagonal 
elements of H are nonnegative and equal to the interaction rates; that is, the 
element Hp with a ^ (3 is equal to the rate of change of the state /3 to the 
state a. a and f3, each represent the state of two adjacent sites. For example if 
a = 01 and (3 = 10, then H% is the rate of particle diffusion to the right. 
Using 

s (g s(a • N) <g (b • N)H = a a bp H af3 lS s <g sN~< <g N s , (12) 

where a and b are arbitrary covectors, one can write down the evolution equa- 
tions of the one- and two-point functions of n^s. It turns out that in the evo- 
lution equation of the one-point function, there are two-point functions, and in 
the evolution-equation of the two-point function, there are three point functions, 
unless the reaction rates satisfy the following constraints [27-29]. 

where 

1 i\a tr«/3 

A- 7 <5 .— S/3 tl M 7< 5 

2 A a lS :=s H^ a l5 . (14) 
It can be seen that one can summarize the constraints (|13|) in the compact form 

iJugiu = Auigu, (15) 

where 

u := ^ , (16) 

and it is obvious that 

s-u = 0. (17) 

Now, consider a system satisfying the constraints (|13|) (or equivalently (|15l0 . 
and take the vector v satisfying 
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E ^]v=0, 

s-v = l, (18) 



, d,e=l 
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and the covector a such that 



au=l, a-v = 0, (19) 

that is, the basis {a, s} is dual to {u, v}. This choice of a makes the evolution 
equation of (a • N) homogeneous. In [28,29], it is shown that the matrix in 
the left-hand side of the first equation in (|18|l . has a left eigenvector with the 
eigenvalue zero. (This left eigenvector is s.) So it does have a right eigenvector 
with the eigenvalue zero as well. That is, there does exist a vector v satisfying 
<|18|) . In fact, one can even find a real vector v satisfying i|18|) . From now on, a 
in © is assumed to satisfy ljl§|) . 

Assume further, that the initial condition is translational-invariant. This 
means that the one-point function is independent of the site, and the two-point 
function depends on only the difference of the sites' numbers. It turns out that 
the evolution equation for the one-point function is 

^ = 0* + ")/, (20) 



where 
and 

Also, taking 
one arrives at 



/ := («<), (21) 
[i =s ® a H u eg) v + a Cg> s H v Cg> u, 

v=s®a.Hv®u + a.®sHu®v. (22) 

Fi := (n k n k+i ), (23) 



^=p(F i ^ 1 +F i+1 )+2vF i , 1 < i < L 
dFi 

-^=pF 2 + {v + X)F l +pf + a, (24) 



where 



p :=a <g>ai/(u®v + v®u), 

a :=acg) aH v Cg> v. (25) 

From the definition 123|l . it is also seen that 

F L+1 - i =F i . (26) 

It is seen that only five parameters enter the evolution equation of the up-to- 

two-point functions, and all of these can be expressed in terms of the matrix 
elements of 

h ■= h + n#n, (27) 
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where II is the permutation matrix. These parameters can be rewritten as 

/U :=s ® aH u ® v 
f :=s ® a/f v (g u 

A :=-a<ga!/u<gu 
p :=a <g a_ff u (g v 

er :=ia<g al/v <g v (28) 

3 Solution of the evolution equations 

The solution to Q2U[) (the evolution equation of the one-point function) is easily 
seen to be 

/(t)=/(0)exp[(/* + i/)t]. (29) 
Putting this in Q24J1. the second equation becomes 

dFi 

= p F 2 + {v + \)Fi + P f(0) exp[0u + v)t] + a. (30) 

This, combined with the first equation of (|24p. and the constraint (|26|l . are 
sufficient to obtain the two-point functions from their initial value. To do so, 
one takes a solution like 

F,(t) = J2 F *E(0)eMEt), (31) 

E 

and puts it in the equations. From the first equation of (|24p. one arrives at 

£; J F lf; (0) = A1 [ J F 1 J _ 1B (0)+F, +1B (0)] + 2^F lf; (0), 1< i < L. (32) 
lj5U|) becomes 

£ Fx B (0) = p F 2 B (0) + {v + X)F 1E (0) + p /(0) + a *„ iB . (33) 

To solve (I32|) . one takes 

F iE (Q) = c B + (34) 

Putting this in J52J}, one arrives at 

F = /i(z + z _1 ) + 2i/. (35) 

The equation for z' is similar, and in fact z' is the inverse of z. Then, using 
(|211), one can write ffi^i as 

F lE (0)=c E (z t + z L+1 - 1 ). (36) 
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Putting this in l|33|) . one can obtain the coefficient eg. 

For £? = p+f,or£J= 0, the equation for ce is a nonhomogeneous one, and 
ce is obtained. For E different to the above values, the equation for is a 
homogeneous one, and only for certain values of E there exist nonzero solutions 
for ce- These values of E are among the eigenvalues of H, of course. Using 
(|36|) . i|33|) . and l|35|l . the condition for ce being nonzero is seen to be 

^ L+1 ^ 2 + z^+D/ 2 ] = (A — v)[z< L -W + z^' 2 }. (37) 

Some of the roots of this equation (for z) are phases (their absolute value is 
one). In the thermodynamic limit [L to infinity) it is easy to find the nonphase 
solutions. It is seen that if a solution has absolute value less than one, then in 
the thermodynamic limit, 

A — v 

and it is obvious that such roots of (|37|1 are real. So, in the thermodynamic limit, 
the energy values E entering the translationally-invariant two point function are 
zero, (/i + v), the values coming from (|35[1 with \z\ = 1, and possibly only one 
other value coming from (|35|l with z satisfying l|38[) . The largest nonzero value 
of E, determines the relaxation time towards the equilibrium. 

One point should be noted. In general, the limit of the largest relaxation 
time of a finite system, as its size tends to infinity, may differ from the relaxation 
time of the infinite system. It can be shown, however, that it is not the case for 
our system. In fact, if one solves the eigenvalue equation for the infinite system, 
one has to omit the periodicity condition l|26l) . and use instead a condition 
that the two-point function does not blow up in the limit that the distance 
between the two points tends to infinity. This means that either z and z' (the 
inverse of z) are unimodular, or in (|34|l there remains only one term, the term 
corresponding to the one with modulus less than one. For the latter case, one 
again recovers l|38|l . So, the spectrum of the infinite system, is in fact equal to 
the limit of the spectrum of the finite system with periodic boundary conditions, 
in the infinite-size limit. It is of course true that if in the initial condition, the 
coefficients of some eigenvectors vanish, then the relaxation time may differ 
from the largest relaxation time. But this happens independent of the size of 
the system. Another case when the relaxation behavior of the infinite system 
differs from the limit that of the finite system, is when the spectrum becomes 
continuous to zero, that is, in the infinite limit system, there is no eigenvalue gap 
between zero and the other part of the spectrum. In this case, the relaxation 
behavior of the infinite system may be a power law, rather than exponential 
decaying. But again this is not the case for the present system. To summarize, 
in the present system the largest relaxation time of the infinite system is equal 
to the limit of of the largest relaxation time of the finite system. 
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4 The parameter space determining the energy- 
spectrum of the two-point functions 

Consider a one-dimensional single-species nearest-neighbor-interacting system, 
for which the evolution equations of up-to-n-point functions are autonomous, 
(we call such systems autonomous.) The Hamiltonian H characterizing such a 
system (hence satisfying Q15[l). contains 10 parameters. As it was seen from the 
previous section, of the parameters entering H, only five parameters enter in 
the evolution equation of the two-point functions. All of these are expressible in 
terms of the symmetrized (with respect to permutation) Hamiltonian H. It is 
easily seen that as H satisfies (|15|) . H satisfies (|15|) as well. So the system char- 
acterized by H, is autonomous as well. Such a system contains 6 independent 
rates. In fact, one can write H as 



H = 



2ri - r 2 




r3 




r-3 


rs 


ri 




- r 7 


- Ti 




r% 


m 








-r3 - rr - r 4 


r% 


r-2 




'4 






-r§ - 




n - 


Vr 2 - 


hr 3 = 


U + r 5 + r e . 





(39) 



with 

(40) 

Of the five parameters entering the evolution equation of the two-point function, 
only three parameters determine the energy-spectrum. These are [i, v, and A: 

/i = r 7 + r 4 - n - r 2 = r 7 + r 3 - r 5 - r 6 , 
v = -r 7 - n - r 2 - r 3 = -r 7 - r 4 - r s - r 6 , 

A = - ri+r3 + r4 + re . (41) 
From these relations, it is seen that 

f<-|A»l<0, 

v < A < 0. (42) 

As the rates are nonnegative, if v = 0, then H = 0, which makes the two-point 
functions constant. Assuming v ^ 0, one can scale time and make v = — 1. 
So, apart from a time-scale, there are only two parameters determining the 
energy-spectrum, fx and A: 

H <i, 

-1 < A < 0, 
v = -1. (43) 

It can be shown that the whole region of the above is physical. That is, corre- 
sponding to any A and fi satisfying the above inequalities, there are autonomous 
systems yielding the desired A and \i. To prove this, first consider four specific 
systems: 
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• 


' - 1 


= r 6 


= 1, 


>2 


= r 3 


= r 4 


= r 5 


= r 7 


= 0, = 


> GM) = 


(-1,-1) 


• 




= r 5 


= 1, 
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= r 3 


= r 4 


= r 6 


= r 7 


= 0, = 


> (H, A) = 


(-1,0). 


• 


rs 


= r 4 


= 1, 


n 


= f2 


= r 5 


= r 6 


= r 7 


= o, = 


> (/i, A) = 


(1,-1). 
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T7 


= 1, 


n = 


T2 


= r 3 


= r 4 


= r 5 


= r 6 


= 0, = 


> (M, A) = 


(1,0). 



Any point (/i, A) in the region described by l|43|l . can be written as 



( M , A) - Cl (-1, -1) + c 2 (-l, 0) + c 3 (l, -1) + c 4 (l, 0), (44) 
where c a 's are nonnegative. A system with the rates 
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n = c a r ai , (45) 

a=l 

where r a i is the rate of the a'th system introduced above, gives the desired 
( M ,A). 

5 Dynamic phase transitions in the two-point 
function 

It was shown in the previous section, that for any nonzero Hamiltonian v ^ 0, 
so that one can normalize v to — 1. In section 3, it was shown that the energies 
entering the two-point function are 0, E\ := fi — 1, and fi(z + z _1 ) — 2, where 
in the thermodynamic limit \z\ — 1, or at most one non-unimodular z exists, 
the value of which comes from (|38|l . This is provided the absolute- value of the 
left-hand side of (|38|l is less than one. So, the energies (apart from 0) are E\, 
any number in the interval Iq := [— 2 — 2\fi\, — 2 + 2|/x|], and possibly 

E 2 :=\-l + ^- l , (46) 

The largest relaxation time of the two-point function is —E~l x , where E max is 
the largest nonzero value of the energy spectrum. The relaxation time of the 
one-point function is —E±. The fact that the energy spectrum of the one-point 
function consists of a single value, is a result of the translational-invariance of 
the initial state of the system. Otherwise, there would be many energies for the 
one-point function, which could lead to a dynamical phase-transition in the one- 
point function [30-33] . The comparison of the relaxation times of the two-point- 
and the one-point-functions, is a comparison of E max and E±. If the former is 
larger, the largest relaxation time of the two-point function is larger than the 
relaxation time of the one-point function (the slow phases). If the two are equal, 
the relaxation-times are equal (the fast phases). So, the relation of E\, E 2 , and 
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Iq , determines the relaxation behavior of the two-point function (its dynamical 
phase). It is seen that 

J < E x , fi> -i, 

Ei £ Jo, M < -\, (47) 

where < 1 has also been used. If E\ > Iq, then the relaxation time of 
the one-point function is equal to the largest relaxation time of the two-point 
function. If E\ £ Io, the the largest relaxation time of the two-point function is 
larger. 

For Ei to be among the energies, the absolute value of the left-hand side of 
should be less than one. So, 

fiE = E 2 , A<H-1, 

3E = E 2 , A > H-l. (48) 

Finally, 

(A + l)(E 2 - E^ = A 2 + /i 2 - A n + A - /i =: f(p, A), (49) 
from which (using A + 1 is nonnegative) , 

E 2 <E 1 , /(a*,A)<0, 

E 2 >E l , f(n,\)>0. (50) 

/ = is an ellipse, the interior points of which correspond to E 2 < E\ , and the 
exterior points of which correspond to E 2 > E\ . 

These three inequalities divide the whole phase space < 1, 1<A<0) 
into five phases: 

I) fi < -i, A < \n\ - 1. 

In this phase, E\ € Iq, and E 2 is not an energy. This is the slower phase, 
and the largest energy is E maK = —2 — 2/i. 

II) M<-|, A>|/x|-l. 

In this phase, E\ £ Iq, and E 2 is an energy, in fact the largest one. This 

2 

is the slowest phase, and the largest energy is E max = — 1 + A + 

III) /i> -|, A < \n\ - 1. 

In this phase, E\ > Iq, and E 2 is not an energy. This is the fastest phase, 
and the largest energy is E max = — 1 + fi. 

iv) M > -|, H - 1< a < ^ 1+ ^ /(1 2 +3M)(1 ^ ) -. 

In this phase, E\ > Iq, E 2 is an energy, and E 2 < E\. This is the fast 
phase, and the largest energy is E max = — 1 + /i. 
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In this phase, E\ > Iq, E% is an energy, and E<i > E\. This is the slow 

2 

phase, and the largest energy is E max = — 1 + A + ^rj. 

This phase structure is summarized in fig. 1. 

As previously mentioned, these phases arise from the fact that the energy- 
spectrum of the two-point function consists of a continuous part, an energy 
equal to the energy appearing in the one-point function, and possibly another 
energy. This shows that the relaxation of the two-point function is at least as 
slow as that of the one-point function, and may be slower, depending on the 
relative position of the discrete and continuous parts of the spectrum. The fast 
phases (phases III and IV), are those in them the relaxation time of the one- 
and two-point functions are equal, while in the slow phases (phases I, II, and 
V), the relaxation of the two-point function is slower than that of the one-point 
function. 



6 A one-parameter family as an example 

Consider a system with the Hamiltonian 
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This Hamiltonian describes a system with the following reactions. 

%A — > any other state, with the rate lo/A, 
A% — > any other state, 
00 — ► any other state, 
AA — > any other state, 
It is seen that for this system, 

H 

and 

2lo, 



with the rate lo/4, 
with the rate (1 — w)/4, 
with the rate (1 — lo)/4. 

2H, 



ft 
A 



(51) 



(52) 
(53) 

(54) 



For this system, p defined through (|28|l is equal to zero, hence there is no term 
proportional to e Elt in the right-hand side of (|30|l . However, it is seen that 
adding a term 

/-I \ 
1 
1 
\l 0-2/ 



Hi 



(55) 
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to the Hamiltonian H in 151|) . changes the values of p, v, fj,, and A to 



P =r, 

v = -l-2r, 

fi = -l-2r + 2u>, 

A = - ~ - r. (56) 

For small values of r, one can use v = 1, and /i and A as in Ij54(l . Then, with 
different values of uj, this system can exist in all the above five phases: 

1 1 

4 <W< 3' 
1 5- V5 

3 8 ' 

5- V5 3 

\<u<l. (57) 

It is seen that increasing to, the system undergoes phase transitions from the 
phase I to II, then V, IV, and finally III. 
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Figure 1: the dynamical phase structure in the (fi, A) plane 
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